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Abstract 

We compute the decay rate of the Standard Model Higgs boson to bottom quarks 
to order aag- We apply the optical theorem and calculate the imaginary part of 
three-loop corrections to the Higgs boson propagator using asymptotic expansions 
in appropriately chosen mass ratios. The corrections of order aag are of the same 
order of magnitude as the 0{ag) QCD corrections but have the opposite sign. 

PACS numbers: 14.80.Bn 14.65.Fy 12.15.Lk 


1 Introduction 


After the discovery of a Higgs boson in rnn I of the CERN Large Hadron Collider it is 
one of the main tasks of run H to determine the properties of the new particle. Among 
them is the coupling to other particles. This is predominantly done by determining Higgs 
production cross sections and decay branching ratios, i.e. the ratio of the partial decay 
width of the Higgs boson to the considered particles normalized to the total decay rate. 
The latter is dominated by the partial decay rate to bottom quark, bb), which 

hence influences all branching ratios. Thus, r{H —)■ bb) should be available as precisely 
as possible. 

QCD corrections are known up to order (see, e.g.. Refs. IDE]) and hrst results of order 
af induced by virtual top quarks have been obtained in Ref. [8]. Good convergence of 
the perturbative series is observed leading to a 0.1% contribution of the corrections 
to r{H —> bb). As far as electroweak corrections are concerned only one-loop corrections 
are available which have been computed beginning of the nineties At two- and 

three-loop order only the leading corrections are available mHH]. In this work we 
compute QCD corrections to the full 0{a) result and thus obtain all contributions of 
order aag to the partial decay rate of a Standard Model Higgs boson into bottom quarks. 
Analog corrections to the decay rates of the Z and W bosons have been computed in 
Refs. [I5] - [T7] and [18], respectively. 

We parametrize the corrections to the decay rate as follows 

r{H bb) = (1 + + ...), (1) 


where the ellipses stand for higher order corrections in a and Og- It is convenient to split 
the electroweak corrections into a weak and a QED contributions which to our order are 
separately hnite and gauge invariant: 

A^“^ = -I- , 

_ ^(QED,Os) _|_ ^(weak.Os) /'2'\ 


The aim of this paper is the computation of the mixed corrections In Eq. ([T]) 

denotes the Born decay rate which is given by 


p(o) 


NcamlMn 3 

8s^M^ 


( 3 ) 


where = 3 is the number of colours and sw is the sine of the weak mixing angle. 
/3o = y/l — 4m ^/Mjj is the velocity of the produced bottom quarks which from now on 
we approximate to /So = 1- As an alternative to Eq. (I3|) one can replace the hue structure 
constant by Fermi’s constant via 


Gf vtq; 1 

\/2 1 - Ar ’ 


( 4 ) 
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where the hnite quantity Ar parametrizes the radiative corrections to the muon decay 
beyond QED corrections within the effective four-fermion theory [19]. This leads to 


p(o) 


NcGpmlMH 

4x/27r 


(5) 


For later reference we provide the Born decay rate including higher order terms in e = 
(4 — d)/2 which are useful in the renormalization procedure. For d ^ 4 both (j3|) and ([5|) 
have to be multiplied by the factoi0 


fie) 



1 + e + e' 



+ 0{e^) 


( 6 ) 


For later convenience we also list the one-loop QCD corrections including terms of order e 
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where Cp = 4/3. To obtain this result the bottom quark mass has been renormalized in 
the MS scheme. The one-loop QFD corrections are obtained from with the help of 

^(QED) ^ ^ 

Cpas 

where Cp = 4/3 and Qh = —1/3 is the charge of the bottom quark. 

The remainder of the paper is organized as follows: In the next Section we discuss the 
method we want to use for the three-loop diagrams of order aag and apply it to the 
one-loop electroweak corrections. The 0{aas) corrections are presented afterwards in 
Section 131 In Section 0] we discuss the numerical effect, compare with the known QCD 
corrections and conclude. 


2 Corrections of order a 

Before discussing the computation of the genuine diagrams of order a we briefly elaborate 
on the counterterm contribution. We follow Ref. [20] and introduce one-loop counterterms 
for the Higgs boson wave function {dZp), the vacuum expectation value {6v) and the 
bottom quark mass (dm/,)- This leads to the following counterterm contribution of 

Ajweak) ^ ^(0) (^1 - 2^ + dZp - Ar + 25^,^ , (9) 

^Throughout this paper we adopt a MS-like convention and set jp and log(47r) to zero. 


3 















where is given in Eq. ([5]). We do not include the on-shell wave function renormal¬ 
ization of the quarks in since it is automatically taken into account when using 

the optical theorem (see below). In the on-shell scheme the mass counterterm is defined 
through 


m» = M, (1 + iSD . (10) 

with Mb being the on-shell mass. We take from Ref. [211 [22] dropping all tadpole 
contributions. The divergent part of determines the MS counterterm 5^® and is in 
our approximation (i.e. at most ml terms in the decay rate) given by 


x(weak),MS 

"mb 
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32slM^ 
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9 -|- 6siy — 8s 






( 11 ) 


The remaining terms on the right-hand side of Eq. ([9]) can be written in the form 
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( 12 ) 


where v is the vacuum expectation value and and denote the two-point 

functions of the corresponding bosons in the notation given in Ref. [20]. The prime in 
the case of the Higgs boson two-point function denotes the derivative w.r.t. the external 
momentum squared, q^. Afterwards = Ml is chosen. 
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Figure 1: Sample Feynman diagram contributing to the 0{a) corrections of r{H bb). 
External dashed lines denote the Higgs boson. 


For the evaluation of the decay rate r{H —)■ bb) we use the optical theorem which for our 
application has the form 

TiH ^ bb) = ^Im [Eniq^ = + *e)] , (13) 

Mr 

where Y^niq^) is the Higgs boson two-point function which is evaluated on the Higgs boson 
mass shell. As a consequence we have to consider Feynman diagrams as shown in Fig. [T] 
to evaluate the 0{a) corrections. In this approach we automatically take into account the 
on-shell wave function renormalization which is the reason why we have not considered it 
in Eq. ([9]). Note that we neglect corrections except the leading ml factor and thus the 
contributing diagrams either contain Z bosons (possibly together with neutral Goldstone 
bosons) or W and/or charged Goldestone bosons and top quarks. 

We express our hnal result in terms of the MS bottom quark mass which, as is well 
known, leads to a better perturbative behaviour of the decay rate. In this context we 
briefly want to discuss the tadpole contributions to the bottom quark propagator (see 
also discussions in Refs. [2T1[22]). In fact, besides the diagrams in Fig. [1] there are also 
contributions where a closed loop is connected via a Z or Higgs boson to the bottom 
quark line, so-called tadpoles. These contributions are exactly canceled by the on-shell 
counterterm contributions to the bottom quark mass. For this reason we drop the tadpoles 
in both parts from the very beginning. Note, however, that after dropping the tadpole 
contribution in the counterterm dm^, it becomes dependent on the electroweak gauge 
parameters ^wjz- The same is true for the contribution from the diagrams in Fig. [T] In 
the sum ^wiz drops out. In case the bottom quark mass is renormalized in the MS scheme 
there is no cancellation and the hnal expression for r{H —> bb) remains .^vv/x-dependent. 
Note, however, that also the numerical value of mi, in the MS scheme (formally) depends 
on ^wjz since in the extraction of mi, from the comparison of theoretical calculations 
and experimental data (see, e.g.. Ref. [23]) no electroweak tadpoles are included. The 
^m/z-dependence in F(Fr ^ 66) and mi, cancels. 

In our calculation we adopt Feynman gauge in the electroweak sector but allow for general 
gauge parameter in gluon propagator. In our hnal result drops out which is a 
welcome check. Our Feynman integrals involve the mass scales Mr, Mt, Mw and Mz- 

Before presenting numerical results let us hx our input parameters which are given by [23F 
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Mt 

= 173.34 GeV, 

Mh 

= 125.09 GeV, 

Mw 

= 80.385 GeV, 

Mz 

= 91.1876 GeV, 

mfe(mfe) 

= 4.163 GeV, 

Gf 

= 1.1663788(6) x 10 

o^s{Mz) 

= 0.1185 


-5 


GeV 


-2 


( 14 ) 


where the four-loop QCD conversion [26] of the on-shell to the MS top quark mass leads 
to mtirrit) = 163.47 GeV and mt{MH) = 166.97 GeV. 

Let us in a hrst step discuss the contribution from the Feynman diagrams which do not 
involve top quarks. As massive particles they only contain Z bosons or neutral Goldstone 
bosons and thus they depend on where = Mjj is the square of the external 

momentum. At 0{a) an exact calculation is possible, however, at 0{aas) the occurring 
integrals become complicated. Thus we evaluate this class of Feynman diagrams in the 
limit M| and apply a Fade approximation to construct an approximation for the 

physical limit = Mfj. In principle one could also imagine to consider q^ ^ M|. How¬ 
ever, this limit contains decays of the form H ZZ which are kinematically forbidden. 
On the other hand, for q^ M| we neglect contributions from H —)■ Zhh, which are, 
however, strongly phase-space suppressed. Furthermore, it is possible to experimentally 
distinguish this hnal state from H —)■ bb. Note that the decay H —)■ Zbb is not included 
in the result of Ref. [10] . 

In the limit q^ <C M| we obtain for the expansion 


(weak,Z) 



(15) 


where the coefficients are given in Table [U includes all relevant contributions 

from and is thus hnite in the limit e —?■ 0. The bottom quark is renormalized in the 

MS scheme. The counterterm contribution is not expanded in and is contained 

in the coefficient Dq. Furthermore, we choose /i^ = Mjj for the renormalization scale. 

In a next step we use the results in Table [T] and construct various Fade approximations. 
The results are shown in Table [2] where also the exact result for from Ref. [10] 

is displayed. The deviation of the numerical approximation based on the [4/4] Fade 
expression and the exact result uni is about 0.01% which justihes the use of this method 
at order aa^. 

Let us next turn to the contribution involving top quarks and a W and/or charged Gold- 
stone bosons. To simplify the integrals and to obtain simple hnal expressions we assume 
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k 

Dk 

0 

-0.008768 

2 

-0.000847 

4 

-0.000196 

6 

+0.000269 

8 

-0.000534 

10 

+0.000894 

12 

-0.001515 

14 

+0.002568 

16 

-0.004383 


Table 1: Coefficients from Eq. ffTSj) for /i^ = Mh- 


Fade approximant 

(weak,Z) 

[3/3] 

-0.009744 

|3/4] 

-0.009746 

|4/3] 

-0.009747 

[4/4] 

-0.009746 

exact 

-0.009747 


Table 2: Numerical results for obtained from the construction of Fade approx¬ 

imations using the coefficients in Table [T] The last row contains the exact result from 

Ref. do]. 


one of the following hierarchies 

(A) Ml < 4M^ < 4M2 , 

(R) Ml < 4M^ ^ 4M2 . (16) 


We stress that the 0(a) corrections can be computed without any assumptions on the 
relative size of the involved mass scales. However, at order aas the hierarchies in Eq. ffT6l) 
signihcantly simplify the calculation. 

In Eqs. (fT6|) the strong hierarchy () means that we apply an asymptotic expansion [28] 
in the corresponding mass ratio. In the case of an approximation sign we Taylor-expand 
the integrand in the mass difference. As a result we obtain in the form 


yy (weak) _ 


i>-l 

E' 

i>-l 


M) 

"21 


SB) 

■'2i 


M2 

M2 




(A 

2i 


i>-\ 


Mw 


M} 




(S) 

2i 


i>-l 


( 17 ) 
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k 

(MS) 


^ (weak) 

-2 

+0.005146 


+0.005146 

0 

-0.004506 - 

0.009746|z 

-0.009106 

2 

-0.000166 


-0.009272 

4 

-0.000100 


-0.009372 

6 

-0.000105 


-0.009477 

8 

-0.000088 


-0.009565 

10 

+0.000048 


-0.009517 

12 

-0.000029 


-0.009546 

14 

+0.000001 


-0.009545 

exact 

-0.009549 

k 

0)^"^ (on-shell) 

^ (weak) 

-2 

+0.004842 


+0.004842 

0 

-0.004754 - 

0.009746U 

-0.009659 

2 

-0.000145 


-0.009804 

4 

-0.000095 


-0.009899 

6 

-0.000078 


-0.009977 

8 

-0.000065 


-0.010042 

10 

+0.000029 


-0.010012 

12 

-0.000018 


-0.010031 

14 

+0.000000 


-0.010030 

exact 

-0.010034 


Table 3: Coefficients C\f^ as defined in Eq. flTTl) (middle colnmn). In the row of the 
right colnmn the snm inclnding the first n terms is shown. On the top part we adopt the 
MS and below the on-shell definition for the top qnark mass. The Z boson contribntion 
is marked by |^. 

where the coefficients and are expansions in M^/(4M^). Note that by definition 
Cq^'^ and contain the contribntion from 

We show onr resnlts for hierarchy (A) in Table |3] adopting again /i^ = Mjj and the MS 
definition for the bottom qnark mass. For the top qnark mass both the MS and on-shell 
mass valne is used. 

Note that ~ l/x as x —?■ 0. For this reason we show in Fig.[2]the quantity 

as a function oi x = and compare the expansion obtained for the hierarchies (A) 

and (B) with the exact result [10] |§ For the plot we use the on-shell definition of the 
top quark mass and set = Mjj. The numerical values are obtained by keeping M\y 
fixed and varying Mf. We take into account expansion terms up to fc = 12 (see Table [3|) 
which corresponds to the expansion depth which is available at order aag (cf. Section [3D. 

^There is a typo in the quantity ATio+n in Eq.(A.2) of [TU]: the minus sign in front of should 
be a plus sign. 













Figure 2: Comparison of as obtained for the hierarchies (A) and (B) with the 

exact result as a function of x = The black (solid) curve shows the exact result, 

the (red) dashed curve the expansion for x —)■ 0 and the (blue) dotted curve the expansion 
around x = 1. The vertical line indicates the experimental result for x ~ 0.215. For the 
renormalization scale of the bottom quark jj? = Mjj has been chosen. Note that 
behaves as log(a;) for a; —)■ 0. 


One observes that for x ^ 0.4 a perfect description is obtained from hierarchy (A) (red, 
dashed curve) and above x ~ 0.4 the result from hierarchy (B) (blue, dotted curve) agrees 
perfectly with the exact result (black line). For the physical value x ~ 0.215 one obtains 

A^eak) _ -0.009525, (18) 

which has to be compared with the results in Table [3] in the lower panel. There is a 
notable deviation of about 5% to the exact result which has its origin in the divergent 
behaviour proportional to 1/x for a; —)■ 0. For this reason we concentrate in Section [3] on 
hierarchy (A). 


3 Corrections of order aag 

In this Section we consider the quantity of Eq. ([T]). An analytic expression for 

A(QED>“s) can easily be obtained from the 0(ag) QCD corrections (see, e.g.. Ref. 
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after adopting the colour factors. It reads 


^(QED,Os) 


Q\ 


2 OlOls 
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691 3 2 
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9 

2C(3) + 




9 . 2 

+ -In^ 
4 



.( 19 ) 


To obtain we proceed as follows: 


• We consider the imaginary part of the three-loop propagator-type diagrams which 
are obtained by dressing the 0{a) diagrams (cf. Fig. [T]for examples) in all possible 
ways with one gluon. This part can be split, in analogy to the 0{a) corrections, into 
a contribution involving Z or Goldstone bosons and into a contribution involving 
W and/or charged Goldstone bosons and top quarks. 

• The bare bottom quark mass in the Born result has to be replaced by the MS renor¬ 
malized counterpart using corrections of order aag- The corresponding counterterm 
is available from Ref. [2^I3U] which we have checked by an independent calculation. 
It is given by 


9mj \ 

1/1 31 27 Smf 

^7 V96 + 7687|; + 256s^ “ 32M^s^ 

Note that contains poles up to order 1/e^ and thus the Born result is 

needed up to order terms. 

• has to be available up to order e and the bottom and top quark masses 

have to be renormalized using one-loop counterterms of 0{as) which are given by 



^(weak,Q!s),MS _ ^ 


■'mi, 


TT^ 


7 


9 


16 128c^ 128s?, 


+ 





m, 


bare 


Af„ 


_i = + i + 

TT \4e 4 Mg 


( 21 ) 


with q = b,t. The corresponding MS counterterm is obtained by dropping the finite 
part on the right-hand side of the above equation. 

• Pas |;o be available up to order e and the bottom quark mass has to be 
renormalized using one-loop counterterms of 0{a) which is given in Eq. (ITT]) . 

• There is a contribution where Vr from Eq. (IT^ multiplies F^^^A*-"®/ Since the latter 
is hnite we do not need the 0{e) part of Vr- On the other hand, since Vr contains 
1/e poles F^^^AOd is needed including 0{e) terms. 

• The fermion-loop contributions to Vr [see Eq. (IT^ j receive two-loop QGD corrections 
which are multiplied by the Born decay rate. Since the fermionic contribution to 
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S^^(O) vanishes only S'^(O) and get correction terms of order aag- We 

compute them in the limit of a heavy top quark and obtain in the MS scheme 


E’^(O) 

~W 




0{aas) 


Ncaas f 
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mfM^ 
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, Mh^ 
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^ L 


19 


72 

503 


24 




+ -T7tC(2) 
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12 


100800 

9523 
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560 


In (A. 
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15876000 ^ 630 
100687 1 


2514758400 
154559 
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2464 
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ifl 
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ln(i^ 

ml 


19423404000 10010 


ln(i^ 

ml 


Mjj 61 
3240 


+ 0 


( M}^ \\ 


where = mt{^). 


( 22 ) 


The individual terms develop poles up to order 1/e^, which cancel in the sum. 

We are now in the position to present results for the order aag corrections. In analogy to 
Eqs. firsj) and flT7|) we introduce 


y^(weak,Q:s,Z) 


^ (weak,as) 



(23) 


where for convenience is added to the coefficient Cl{^\ 

In the case of we proceed as at order a: we compute nine expansion terms for 

the (formal) limit -C M| and set = Mjj. After including the corresponding coun¬ 
terterm contributions we obtain the expansion coefficients listed in Table 01 Afterwards 
we construct several Fade approximants and obtain the results in Table [5l We observe a 
similar stability as at 0{a) and estimate the hnal result as 


= -0.00195(1), 


(24) 


which has an uncertainty of about 0.5%, an accuracy sufficient for all foreseeable appli¬ 
cations. 

In Table [6] we present the results for the coefficients We observe a continuous 

decrease of the magnitude leading to a stable result with two significant digits after 
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k 

Dk 

0 

-0.001692 

2 

-0.000032 

4 

-0.000547 

6 

+0.000836 

8 

-0.001347 

10 

+0.002183 

12 

-0.003604 

14 

+0.006033 

16 

-0.010223 


Table 4: Coefficients from Eq. fl2^ for /i^ = Mjj. 


Fade approximant 

^(weak,as,Z) 

[3/3] 

-0.001955 

13/41 

-0.001954 

[4/3] 

-0.001960 

[4/4] 

-0.001953 


Table 5: Numerical results for (weak,as,z) obtained from the construction of Fade approx¬ 
imations using the coefficients in Table IH 


including six expansion terms, a similar behaviour as at order a. The seventh and eighth 
terms confirm this approximation. It is also interesting to note that the contribution from 
the Z boson diagrams amounts to about 65% of the total result. Furthermore, the leading 
contribution amounts to less than 20% of but to more than 50% of the W 

boson diagrams, i.e., 


4 Numerical results and conclusions 

In Table [7] we summarize our results for the 0{a) and 0{aas) corrections where the 
electroweak part is split into QED and weak corrections. The contribution from the Z 
boson diagrams is listed for completeness; their contribution is contained in A*^"®®"^^ and 
A("®®k,"s)^ respectively. For comparison also the QCD corrections up to 0(ag) [6] based 
on computations of the imaginary part of the massless Higgs correlators are shown in 
Table [71 Top quark induced QCD corrections due to an effective Hbb coupling, which are 
in general small (see, e.g. Eq. (14) of Ref. [S]), are not shown. 

Both at one- and two-loop order the weak corrections are negative whereas the QED 
corrections are positive. One furthermore observes that the weak corrections are about 
an order of magnitude larger than the QED terms. For = Mjj the weak corrections 
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k 

c'r' (MS) 


^ (weak,Os) 

-2 

-0.000479 


-0.000479 

0 

-0.000514 - 

0.001953|z 

-0.002946 

2 

-0.000044 


-0.002990 

4 

+0.000018 


-0.002972 

6 

+0.000003 


-0.002970 

8 

+0.000005 


-0.002964 

10 

+0.000004 


-0.002960 

12 

+0.000002 


-0.002959 

k 

(on-shell) 

^ (weak,Os) 

-2 

-0.000481 


-0.000481 

0 

-0.000382 - 

0.001953U 

-0.002816 

2 

-0.000032 


-0.002848 

4 

-0.000006 


-0.002854 

6 

-0.000008 


-0.002862 

8 

+0.000011 


-0.002851 

10 

-0.000010 


-0.002861 

12 

+0.000002 


-0.002860 


Table 6: Coefficients at order aag as defined in Eq. 0231) . In the row of the right 
column the sum including the first n terms is shown. On the top part we adopt the MS 
and below the on-shell definition for the top quark mass. 


QCD 

A("=) 

0.2040 

A("s) 

0.0378 

A(“b) 

0.0020 

A("s) 

-0.0014 

QED/QCD 

A(Qed) 

0.0011 

^(QED,as) 

0.0001 




^ (weak) 

^ (weak, Os) 

^ (weak,Z) 

^ (weak, Os, Z) 

weak/QCD 

-0.0100 

-0.0029 

-0.0097 

-0.0020 


Table 7: Numerical result for the QCD, QED and weak one-loop and mixed two-loop 
corrections for = Mjj. Note that and contain the contributions from 

A(we^k,z) ^(weak,as,z)^ respectively. 


amount to about —1% which is significantly smaller than the one-loop QCD correction 
(-1-20%), however, it is of the same order of magnitude as the two-loop corrections 

obtained from the massless Higgs correlator, which amount to -|-3.8% (see, e.g.. Ref. i)- 
At the same value of p the correction term amounts to about —0.3% which is a 

factor three larger than the one-loop QED corrections and which is of the same order of 
magnitude, but with the opposite sign, as the three-loop QCD corrections. It is interesting 
to note that the four-loop QCD corrections are —0.1%. 
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Finally, it is interesting to comment on the assumption the QED and QCD corrections 
factorize, an approach often chosen in case 0{aas) terms are missing. To do this we 
dehne 


^(aos,non-fact.) _ ^(oos) _ ^(o)^(os) 


(25) 


which shall be small in case the factorization approach works. From the numbers in 
Table [7] we obtain 


non-fact.) ^ -0.000831, (26) 

which corresponds to about 30% of 

To summarize, in this letter we have computed the complete 0{aas) mixed corections to 
the decay rate bb). They provide a negative shift of about —0.3% to —)■ bb) 

which corresponds to about 30% of the one-loop electroweak corrections and which is of 
the same order of magnitude as the three-loop QCD corrections. 
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